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ABSTRACT 



The ability to resolve all processes which drive galaxy formation is one of the 
most fundamental goals in extragalactic astronomy. While star formation rates and 
r \ ' the merger history are now being measured with increasingly high certainty, the role of 

gas accretion from the intergalactic medium in triggering star formation still remains 
largely unknown. We present in this paper indirect evidence for the accretion of gas into 
massive galaxies with M* > 10 11 M Q at redshifts 1.5 < z < 3 using results from the 
q ! GOODS NICMOS Survey (GNS). Our method utilises the observed star formation 

$— i 1 rates of these massive galaxies based on UV and far-infrared observations, and the 

I amount of stellar and gas mass added due to observed major and minor mergers to 

calculate the evolution of stellar mass in these systems. We show that the measured gas 
mass fractions are inconsistent with the observed star formation history for the same 
galaxy population. We further demonstrate that this additional gas mass cannot be 
■ accounted for by cold gas delivered through minor and major mergers. We argue that 

V/~j | to sustain star formation at the observed rates there must be additional methods for 

Q\ . increasing the cold gas mass, and that the likeliest method for establishing this supply 

Q\ • of gas is by accretion from the intergalactic medium. We calculate that the average 

\ gas mass accretion rate into these massive galaxies, which is later turned into stars 

\Q ■ between 1.5 < z < 3.0, is M — 83 ± 36 M Q yr" 1 . This is similar to what is predicted 

f^*) ' in detailed simulations of galaxy formation. We show that during this epoch, and for 

C\J . these very massive galaxies, 61±21% of stellar assembly is a result of gas accretion, 

while the remaining ~ 39% is put into place through mergers. This reveals that for 
the most massive galaxies at 1.5 < z < 3 gas accretion is the dominant method for 
instigating galaxy formation. 



INTRODUCTION galaxies must form in a hierarchical way through mergers 



1 and accretion of material from the intergalactic medium. 

■ Both observations and theoretical models now overwhelm- Thig vigw hag supporting evidence from both observations 

ingly suggest that galaxies in the universe have evolved from of the magg functions of galaxies (e g ^ Bundy et al 

an early galaxy population dominated by lower mass sys- 2006; Mortlock et aL 20 11), direct observations of mergers 

terns undergoing significant star formation in the early uni- fa thfi digtfmt universe (e g ^ Conselice et aL 20 03, 2008; 

verse, to the large and relatively passive galaxies that we Bluck et &J 20Q9] 2Q12)] ag weU ag through the evolution of 

find today (e.g., Conselice 2006; Bouwens et al. 2010). How galaxy gizeg (e g Ferguson et aL 2004 . TmjiUo et aL 2007; 

this transformation occurs, that is how we get from young Buitrago et al. 2008) . Evidence from internal kinematics also 

low mass galaxies to the large massive galaxies we see in suggest s that galaxy formation is driven at least in part 

today's universe, is a highly debated topic. Essentially, we by the accumulation of gas from the intergalactic medium 

want to answer the question - how do galaxies assemble their faUing Qnto a galaxy (e g Kereg et al 2005; Qenzel et aL 

stellar masses? The answer to this question will have pro- 200g; Dekd et ^ 200g . Bournaud et al 2011) 
found implications for both the physics of galaxy formation 

and for understanding properties of the universe itself. While it is now generally accepted that galaxy forma- 

Historically, it was once thought that galaxies formed in tion is a hierarchical process driven by the accumulation 

a manner similar to stars through a collapse of gas that later, of stars and gas located outside of a galaxy after its ini- 

through some process, underwent intense star formation. In tial formation, the details of this assembly remain largely 

this paradigm the total baryonic mass of a galaxy does not unknown. Cold gas infall into galaxies, otherwise known as 

change significantly with time, and its stellar mass evolves 'cold gas accretion' is theorised to be an important aspect in 

by rapidly converting gas into stars. However, it is clear that the process of massive galaxy formation (e.g., White & Frenk 
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1991; Birnboim & Dekel 2003; Keres et al. 2005, Dekel et al. 
2009), and even perhaps the dominant method by which 
galaxies assemble their mass. However, there is little to no 
evidence for gas accretion onto galaxies at high redshift cur- 
rently (e.g., Steidel et al. 2010), although some claims are 
appearing at lower redshifts (e.g., Rauch et al. 2011). This 
is primarily due to the difficultly of observing this process 
since the covering fraction of accreting cold gas is likely not 
high (e.g., Dekel et al. 2009), nor would this gas be easily de- 
tected in, for example, absorption (e.g., Weiner et al. 2009; 
Giavalisco et al. 2011). 

In this paper we develop a new approach to this problem 
by analyzing the evolution of the massive galaxy population 
in terms of the history of its baryonic assembly. We exam- 
ine how the stellar mass of a galaxy is built up over time 
through various galaxy formation processes that can be ob- 
served within a unique sample of M, > 10 11 Mq galaxies at 
1.5 < z < 3 taken from the GOODS NICMOS Survey (GNS) 
(Conselice et al. 2011). By examining the addition of stellar 
mass due to major and minor mergers (Bluck et al. 2009; 
2012), and the observed star formation history (Bauer et al. 
2011) we provide through this method circumstantial evi- 
dence for gas inflow, or gas imported in through extremely 
minor mergers, as an important process in galaxy assembly. 

After comparing the amount of accreted mass to the 
mass assembled through merging, we conclude that gas ac- 
cretion is likely the dominant process for massive galaxy 
assembly at this epoch. We further discuss the comparison 
with theory, and describe some of the implications of our 
results, including how our work relates to the G-dwarf prob- 
lem (e.g., Larson 1974) the rapid gas depletion time-scales 
in galaxies. 

This paper is organised as follows: §2 includes a discus- 
sion of the data sources we use in this paper, and the sample 
selection, §3 is a description of our baryonic mass assembly 
analysis, §4 presents our arguments for gas accretion and §5 
is our summary. We use a standard cosmology of Ho = 70 
km s" 1 Mpc -1 , and f2 m = 1 — fl\ = 0.3 throughout. 



2 DATA 

The data and methods we use in this paper originate from 
the GOODS NICMOS Survey (GNS), and are described in 
detail in Conselice et al. (2011), Mortlock et al. (2011), Bluck 
et al. (2009, 2012) and Bauer et al. (2011). The primary 
galaxies we examine are 81 massive galaxies with stellar 
masses M* > 10 11 M at redshifts 1.5 < z < 3. These were 
selected through a variety of colour selections, include Dis- 
tant Red Galaxies (DRGs), IRAC Extremely Red Objects 
(IEROs) and the BzK galaxies (Conselice et al. 2011). We 
also utilise photometric redshifts which are described in de- 
tail in Griitzbauch et al. (2011a,b). 

To calculate stellar masses we use a Bayesian method 
to fit spectral energy distributions based on various star for- 
mation models to the galaxy photometry. The distribution 
of the resulting stellar masses gives us an error in their mea- 
surement, and we take the peak value of the distribution 
as a measure of our stellar mass. A Salpeter IMF was used 
in these calculations. Further details about our sample se- 
lected are included in Conselice et al. (2011) where we refer 
interested readers for more details. 
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Figure 1. The relationship between the star formation rate and 
effective radius, R c , divided into three redshift ranges between 
1.5 < z < 3 for galaxies with stellar masses M, > 10 11 Mq . 



The galaxy effective radii which we use in this study to 
measure star formation rate densities, originate from GAL- 
FIT fits to our NICMOS imaging, taken from the methods 
and catalog of Buitrago et al. (2008, 2011). These are single 
Sersic profile fits to the light profile in which we retrieve both 
a Sersic index, n, as well as the effective radius, R c . These 
sizes are measured solely using the GOODS NICMOS Sur- 
vey Hi6o-band imaging. Detailed simulations show that we 
are able to retrieve these parameters given the observing 
conditions of our sample (Buitrago et al. 2011). 

We also utilise results from previously published GNS 
studies which describe the evolution of this sample in terms 
of the star formation rate (Bauer et al. 2011) and the merger 
history (Bluck et al. 2009; 2012) as a function of stellar mass. 
We give brief summaries of these analyses in reference to our 
current study throughout this paper. The star formation 
rates we use for our galaxies are measured through rest- 
frame UV and Spitzer 24 /jm measurements. For the UV 
star formation rates which we use for most of the analysis 
here, we utilise ACS z-band imaging to obtain a rest-frame 
UV flux measurement which is converted to a star formation 
rate after applying a k-correction (Bauer et al. 2011). We 
further measure the UV slope for each galaxy using the UV 
colour to correct for dust extinction. 

Our galaxy sample is chosen simply as all systems 
within our sample which have stellar masses of M* > 10 11 
Mq , and are at redshifts 1.5 < z < 3. We do not select 
certain types of galaxies within this selection, such as pas- 
sive galaxies or star forming galaxies. Although amongst this 
sample less than 5% are passive galaxies (Bauer et al. 2011), 
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we include all of the galaxies in this selection as the star 
formation rate and gas mass fractions we use are an aver- 
age within our stellar mass selection. We assume that our 
sample of galaxies represents different phases of a massive 
galaxy's life-span throughout the redshifts 1.5 < z < 3 - 
for example, the fraction of galaxies which are passive rep- 
resents the fraction of time a typical galaxy in our sample 
is not undergoing star formation. Our sample of galaxies all 
have very similar stellar masses and environments, suggest- 
ing that these systems should have very similar formation 
histories (e.g., Griitzbauch et al. 2011a,b). The result of this 
is that we can obtain the average amount of stellar mass 
added to this sample selection due to star formation and 
merging without having to worry about biases from differ- 
ent evolutionary phases. 

This selection also limits our ability to trace the exact 
same galaxies over redshifts due to the fact that galaxies are 
growing in stellar mass due to merging and star formation, 
and thus there will be more galaxies in our stellar mass limit 
towards lower parts of the redshift selection. In Conselice et 
al. (2011) and Mortlock et al. (2011) we discuss the stellar 
mass evolution for our sample. By examining our mass se- 
lection used in this paper we find that the number density 
of galaxies selected by the criterion M* > 10 11 Mq grows 
by a factor of between two and four from redshifts z = 3 to 
2 = 1.5. We examine later in this paper how our selection 
would change if we examine a massive galaxy sample at a 
constant co-moving volume (e.g., Papovich et al. 2011). 

Figure 1 shows a summary plot of our data and massive 
galaxy sample. Plotted are the star formation rates vs. the 
effective radii for each galaxy in our sample, divided into dif- 
ferent redshift bins. The average star formation rate for our 
sample is roughly constant with redshift (Bauer et al. 2011), 
and there is a slight evolution in the sizes of these galaxies, 
such that they are growing with time (e.g., Buitrago et al. 
2008). 



3 EVOLUTION OF GALAXY MASS 

3.1 Galaxy Gaseous Mass at z > 1.5 

The evolution of the cold gas mass within galaxies is an 
important aspect for understanding galaxy evolution, yet 
our observations of this quantity are still very basic and are 
prone to systematic errors. Furthermore, all of these mea- 
surements are indirect, even when using CO as a proxy. We 
present in this section arguments for how to quantify gas 
masses in our galaxy sample based on the inverse Schmidt- 
Kennicutt relation (e.g., Kennicutt 1998; Erb et al. 2006; 
Daddi et al. 2010; Boquien et al. 2011) using our measured 
SFRs and sizes (Bauer et al. 2011, Buitrago et al. 2011; §2). 

Obtaining a measure of the cold gas mass in galaxies 
is challenging, and there are several attempts to measure 
this at high redshift (e.g., Erb et al. 2006; Mannucci et al. 
2009; Daddi et al. 2010), yet these methods contain several 
important limitations which must be considered. These in- 
ferred cold gas mass fractions are measured through the in- 
verse Schmidt-Kennicutt relation, which relates the cold gas 
mass density to the star formation rate density (e.g., Ken- 
nicutt 1998). This relation, for nearby galaxies, has been 
studied in detail in recent papers such as Bigiel et al. (2011) 




log (M.) 



Figure 2. The relationship between stellar mass (M„) and the gas 
mass fraction ^t ga s , which is the ratio of the cold gas mass divided 
by the stellar+gas mass. The red crosses are for z > 2 galaxies 
from Erb et al. (2006) and the blue circles are from Mannucci et 
al. (2009). The solid boxes are from the GNS (this work). The 
error bars show the averages and the 1 a dispersion as a function 
of stellar mass. 



who find a significant scatter, but still a correlation, between 
these two quantities using molecular gas and star formation 
measurements. Much of this scatter is however proposed by 
Krumholz et al. (2012) to be due to observational projec- 
tion effects, and that there is a universal star formation-gas 
density law that applies at all redshifts (see also Narayanan 
et al. 2011). There is also evidence that at least two star 
formation laws may apply at higher redshifts, one for disk- 
like galaxies, and another for starbursts (e.g., Bouche et al. 
2007; Daddi et al. 2010; Genzel et al. 2010). 

Most high-z observations, including ones based on CO 
detections, argue that the cold gas mass fraction, M g /M* is 
between 10-50% for galaxies more massive than M* = 10 10 
M Q (e.g., Genzel et al. 2010; Daddi et al. 2010; Figure 2). 
This cold gas mass fraction tends to rise for lower mass 
galaxies (e.g., Erb et al. 2006; Mannucci et al. 2009). There 
is also evidence that the efficiency of star formation within 
starbursts, such as ULIRGs and sub-mm galaxies is more ef- 
ficient than that given by the standard Schmidt-Kennicutt 
law, resulting in lower derived gas mass fractions at a mea- 
sured star formation rate density (e.g., Daddi et al. 2010; 
Boquien et al. 2011). 

We calculate gas masses for the massive M* > 
10 11 Mq sample, and derive the cold gas mass fraction using 
a form of the global Schmidt-Kennicutt law calibrated for 
nearby star forming galaxies. The relation we use is: 

/ y xl.4±0.15 

Esfr = (2.5±0.7) x 10~ 4 gas _ 2 M yr^kpc" 
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Figure 3. The relationship between the gas mass fraction and 
redshift for galaxies with stellar masses M* > 10 Mq . The error 
bars show the 1 a scatter for the average values which are plotted 
as large solid triangles. The red triangles are for this method when 
examining galaxies at a stellar mass cut off of log M* > 11, and 
the open blue triangles are when holding at a constant co-moving 
volume to account for more galaxies entering our mass selection 
over time. 



Figure 4. The gas depletion time scale shown as a function of 
redshift for each galaxy within our sample. The time-scale is in 
Gyr, and the triangles show the average values at different redshift 
ranges, and the errors bars are the la scatter. 



where Esfr is the surface density of star formation, and 
E gas is the surface density of cold gas. We assume in this 
calculation that the star formation follows the distribution 
of H-band light, which is found to be the case (Ownsworth 
et al. 2012 in prep). We calculate the star formation rate 
density within each galaxy based on the effective radius, R e 
(§2), and half of the measured total star formation rate. 
From this we obtain the gas mass density using eq. (1) and 
we then calculate the total gas masses in our systems based 
on this. 

In Figure 2 we plot the gas mass divided by total bary- 
onic mass (M* + M g ) for our sample massive galaxies. We 
also show in Figure 2 comparison values from Erb et al. 
(2006) who study galaxies at similar redshifts, but at lower 
masses, and Mannucci et al (2009) who study similar star 
forming galaxies at z > 2. Figure 2 shows that the galaxies 
in our sample, as shown by the black boxes, have gas total 
mass fractions which are around 10%, with some scatter. 
However, our values are lower than the gas total mass frac- 
tion found for lower stellar mass systems, and extend the 
trend found in Erb et al. (2006) and Mannucci et al. (2009), 
such that the highest mass systems have the lowest relative 
gas fractions. 

To determine how much more stellar mass these galax- 
ies could acquire from the existing cold gas in star formation 
events, we show in Figure 3 the gas mass fraction defined 
as the cold gas mass divided by the total stellar mass. The 
averages at the three redshift ranges and their la disper- 
sions are at 2 ~ 1.75 - f s = 0.13, a = 0.15, z ~ 2.29 - 



/ g = 0.08, a = 0.08 and at 



2.68 - f e = 0.17, a = 0.19. 



We also show these averages if we hold our selection at a 
constant co-moving volume as the open blue triangles. We 
thus find, as others have, that the gas mass density is roughly 
constant with redshift, as is also found by comparing deriva- 
tions of gas mass fractions at different redshifts for lower 
mass galaxies (e.g., Erb et al. 2006; Mannucci et al. 2009). 
Furthermore, these average values not only have a large dis- 
persion, but also contain potential systematics which we ex- 
amine later. We find very similar gas mass fractions when 
we follow a constant number density of galaxies from high to 
low redshift within 1.5 < z < 3. This is the result of the gas 
mass fractions being largely similar when using a M* > 10 11 
Mq selection, or when using a constant number density se- 
lection, as used in e.g., Papovich et al. (2011). This is due to 
the fact that the constant number density selection is always 
within the M* > 10 11 Mq selection, and within this mass se- 
lection the values of the gas mass fraction are uniform on 
average. 

It has been noted in previous studies that the gas masses 
observed are not high enough to account for the prolonged 
star formation in these systems (e.g., Genzel et al. 2009). 
Figure 4 displays the depletion time-scale, the gas mass di- 
vided by the star formation rate (ip), i.e., M. s /ip - showing 
that nearly all galaxies in our sample would be depleted of 
gas in less than 0.5 Gyr assuming 100% efficiency of star for- 
mation. This is much less than the time within the redshift 
interval we examine in this paper. This suggests that some- 
thing is replenishing the gas within these massive galaxies 
over time. 
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3.2 The Evolution of Stellar Mass 

Galaxy stellar mass is a measure of both the amount of 
gas which has been converted into stars over time, as well 
as the amount of mass which has been accreted into the 
galaxy from previously existing galaxies. We are able to trace 
both of these properties for the most massive galaxies in the 
universe up to z ~ 3 (e.g., Conselice et al. 2011; Mortlock 
et al. 2011; Bauer et al. 2011; Bluck et al. 2012). 

The reason we can do this now is that we have con- 
structed a complete sample of massive galaxies at redshifts 
2 = 1.5 to 2 = 3 which have stellar masses in excess of 
M» = 10 11 Mq (Conselice et al. 2011). By examining galax- 
ies at such a high mass limit we are ensuring to some de- 
gree that we are examining similar galaxies at different red- 
shifts, given that very few galaxies have masses much higher 
than this. Although some galaxies just below the M* = 10 11 
Mq border at the high end of our redshift range will spill 
over into it during the interval 1.5 < z < 3, the properties 
of these slightly lower mass systems are similar to those at 
M» > 10 11 M (Conselice et al. 2011). Furthermore, all of 
the quantities we discuss in this paper are averages, rather 
than integrated or total properties, which limits the effect 
of having new galaxies appear in the selection. 

In previous work, we carried out detailed calculations 
of the merger and star formation histories for a sample of 
81 massive galaxies at 1.5 < z < 3 (Bauer et al. 2011; Bluck 
et al. 2012) which we discuss here. We find that the average 
star formation rate for galaxies with stellar mass, M» > 
10 11 M is < ip >= 103±8 M yr _1 (Bauer et al. 2011), 
with a dispersion of 75 M yr _1 . We find very similar star 
formation rates for the same sample when we examine their 
Herschel far-infrared spectral energy distributions (Hilton et 
al. 2012 submitted). 

The observations we use in this paper are measured 
from redshifts z = 3toz = 1.5. This time period if roughly 
2.16 Gyrs. We calculate how much stellar mass is created 
during this redshift interval by simply integrating the star 
formation rate. When we do this, we obtain the average 
amount of stellar mass created per galaxy over this redshift 
interval, the value of which is: < %j) > St ~ (2.2 ±0.2) x 10 11 
M . This is roughly the average amount of stellar mass 
originally in our sample of galaxies at our highest redshift, 
z ~ 3, which we denote as M*(t=0) = M»(0). Thus we can 
measure the ratio of the stellar mass created in star forma- 
tion within these galaxies divided by the original stellar mass 
within these systems, or < ip > <5i/M»(0) = 1.1 ± 0.13. We 
do not consider directly here the amount of matter brought 
back into the interstellar medium of these galaxies due to SN 
and other stellar evolutionary events. The amount of this is 
less than 5% of the stellar mass created (e.g., Papovich et 
al. 2011), and is much lower than our systematic errors. In 
any case, it is the possibility that this gas is reused in star 
formation that has implications for our results, rather than 
the removal of stellar mass. 

When considering the accretion of stellar mass into our 
sample of galaxies, we also have to consider both major and 
minor mergers, as another important route for obtaining 
stellar mass. The amount of stellar mass added to a galaxy 
due to the merger process is given by the integral over the 
merger history, based on the fraction of galaxies merging, 
and the likely time-scale for mergers (e.g., Conselice et al. 



2003; Conselice 2009; Bluck et al. 2009; 2012). We carry out 
this integration using the observed merger history and time- 
scale for mergers (e.g., Bluck et al. 2009; 2012). When we do 
this calculation over our integral of time between z = 3 and 
z — 1.5 (see §4.3) we find that the total amount of stellar 
mass added by the merger process is, as a ratio of the initial 
stellar mass, M,,m(£)/M*(0) = 0.54 ± 0.23. 

This implies that the merger process, through both mi- 
nor and major mergers, increase the stellar mass of selected 
galaxies by ~ 50% during the interval z — 1.5 — 3 (see Bluck 
et al. 2012 for details). By making assumptions about the 
gas content of these galaxies, based on the results in §3.1 
we can then calculate how much gas is added over time due 
to these mergers. Furthermore, we can test whether the in- 
ferred amount of gas brought in with these mergers is enough 
to fuel the star formation we see in these systems during the 
epoch 1.5 < z < 3. 



4 COLD GAS ACCRETION 
4.1 Formulation of Problem 

A galaxy has several mass components that interlink with 
each other. These are the stellar mass (M„ ), gas mass (M 9 ) 
and dark matter mass (Mbm ). For the purposes of this pa- 
per we consider that the gas mass is cold. Although hot gas 
is present within galaxies at least in the local universe, the 
cooling time for these galaxies is very long and negligible for 
galaxies with masses similar to the ones we examine in this 
paper. 

The total mass (M to t) of a galaxy can then be written 

as: 

M tot = M» + M g + M DM (2) 

where these terms are all expressible as a function of time. 
We first calculate using these terms the amount of stellar 
mass within a galaxy as a function of time M,(t). As the 
stellar mass in a galaxy will not significantly decrease, the 
value of M* can in general only increase. This increase can 
occur through the accretion of outside material from merg- 
ers (M* t M(t)), and through star formation. We write the 
increase due to star formation as < ip > St, where < ip > is 
the average star formation rate within the galaxy between 
an amount of time St (§3.2). 

We then write the change of stellar mass as a function 
of time, M*(t), as 

M„(t) = M*(0) + M„, M (t)+ < i> > St, (3) 

where M»(0) is the initial galaxy stellar mass, and St = t— to. 
We can further express how the cold gas mass (M g ) changes 
in a similar way, although the cold gas mass can increase or 
decrease with time depending on the star formation and gas 
accretion rate. 

M g (i) = M g (0) + M g , M (t) + M glA (t)- <ip> St (4) 

where in this case we have a term for the amount of cold 
gas brought in from mergers, M g ,M(£), the amount from ac- 
cretion of gas without other galaxies, or cold gas accretion, 

Note that this this is not strictly true, as this gas in principle 
could arise from extreme minor mergers, those with mass ratios 



© 2002 RAS, MNRAS 000, [TH9] 



6 Conselice et al. 



M g; A(i)- Furthermore the star formation ongoing within this 
fiducial system will lower the amount of cold gas over time 
by the amount — < ip > St. We furthermore do not consider 
in these equations that some stellar mass created is later lost 
through stellar evolution processes in SN, nor do we consider 
the loss of gas through outflows which can be significant at 
these redshifts (e.g., Steidel et al. 2010; Faucher-Giguere et 
al. 2011). However, both of these effects would act to ulti- 
mately only increase the amount of gas accretion necessary. 

One of the observations of our sample which we de- 
scribed in §3.1 is that the gas mass fraction, M g /M* is 
roughly constant over the interval 1.5 < z < 3 (e.g., Erb 
et al. 2006; Genzel et al. 2009). In this case, we can write 



(5) 



M g (t) M E (0) 
M„(t) ~ M,(0)' 

which is true in the redshift range we examine. Note that 
by using this condition the actual value of the gas to stellar 
mass ratios is immaterial for the following derivation. Thus, 
if we use a different relation between star formation density 
and gas density we would obtain the same results for the ba- 
sic equations, although the final numerical result may differ 
(§4.3). We later examine how our net results would change 
if there was an evolution such that the gas mass density 
declined at lower redshifts (§4.3). 

Using this measurement, we then write the amount of 
gas mass at a later time t, M g (i) by using equation 5. We 
then expressing the stellar mass evolution with time by using 
equation 3 to replace M*(t) from equation 5 to obtain 



M g (t) = M g (0) x 



M», M (t) <ip> St 



M,(0) 



+ 



M*(0) 



(6) 



where the gas mass evolution is expressed in terms of the 
initial gas mass and as a function of how much relative mass 
is added to a galaxy over time due to the increase from 
star formation and the increase from mergers. These are 
much easier to determine than the corresponding absolute 
amounts (§3.2). 

Therefore using equation (6), and the results presented 
in §3.2, we can write the total gas mass after a time St with 
respect to the initial gas mass, M 9 (0) as 



M g (t) ~ (2.64 ± 0.25) x M g (0), 



(7) 



such that the total gas mass within our galaxies at z ~ 1.5 
must be on average ~2.64 times larger than the gas mass at 
the start of our redshift epoch near z ~ 3 due to the increase 
in stellar mass. 

If we then equate this to equation (4) we can calculate 
that the accreted (or extra) gas mass is given by: 

M g , A (i) = (1-64 ± 0.21) x M g (0)+ < ip > St - M gM (t). (8) 

if we divide this equation by the average initial stellar mass 
M*(0) we get: 

M glA (t) _ (1-64 ± 0.21) x M g (0) , < ip > St M gM (t) 



M*(0) 



M»(Q) 



M»(0) M*(0) 



which are less to 1:100 and which we are not sensitive to in our 
observations (Bluck et al. 2012). However, the total amount of gas 
needed implies that these extreme minor mergers would vastly 
outnumber the visible minor mergers at a level much higher than 
any reasonable extrapolation. 



Therefore by knowing the initial average gas mass fraction, 
the fraction of stellar mass increase in the form of stars, and 
the gas mass fraction brought in through merging, we can 
determine the amount of gas mass, relative to our average 
galaxy's initial stellar mass, brought in through 'pure' accre- 
tion events, i.e., not gas brought in with merging galaxies. 



4.2 Gas Accretion Mass Fraction 

As we discuss in §3.1 we find that the initial gas mass frac- 
tion at z = 3 is f g = 0.17, a = 0.19, and that the frac- 
tion of stellar mass formed in star formation is < ip > 
<5t/M*(0) = 1.1 ± 0.13, and the amount added from merg- 
ers is M*, M /M,(0) = 0.54 ±0.23. The most difficult of these 
values to determine is the ratio of the cold gas mass brought 
in from mergers to the total amount of initial stellar mass 
(M gl M(i)/M,(0)). However, we are able to make some mea- 
surement of this based on the amount of mass accreted from 
mergers (M gj M) and the mass weighted gas mass fraction 
for these systems f g . This gas contribution from the minor 
mergers is our key new measurement that allows us to later 
constrain how much gas is coming from the intergalactic 
medium to form galaxies. 

To do this, we first examine the amount of stellar mass 
added due to the major and minor merger process. As de- 
tailed in Bluck et al. (2012) the merger fraction can be pa- 
rameterized as a function of both stellar mass and redshift. 
Bluck et al. (2012) find that at 2.3 < z < 3 the merger frac- 
tion dependence on mass stellar for galaxies more massive 
than M» > 10 11 M Q is given by: 



/ m (M») = (0.28+/-0.17) x 51og(M t ) a91+/ -°- 35 



(10) 



and at 1.7 < z < 2.3 we find a merger fraction dependence 
on mass range for massive galaxies of: 



f m (M,) = (0.12+/-0.09) x 51og(M*) 



1.24+/-0.37 



(11) 



The total amount of stellar mass accreted is then a double 
integral over the redshift range of interest and over the stel- 
lar masses which we probe, which for the GNS is sensitive 
down to M» = 10 9 - 5 M can be written as, 



Mt.M = 



M 2 



(12) 



Where r m is the merger time-scale, which likely depends on 
the stellar mass of the merging pair (Bluck et al. 2012). 
Our integration as mentioned in §3.2, gives a value of 
M*, M /M.(0) = 0.54 ± 0.23. We integrate the amount of gas 
mass added in a similar way, through understanding what 
the ratio of stellar to gas mass is based on the data shown in 
Figure 2. This allows us to compute what the total amount of 
gas and the total amount of stellar mass accreted from these 
minor mergers is. We show in Figure 5 the relative amounts 
of gaseous and stellar mass as a function of the stellar mass 
of the merging galaxy. This shows that the contribution of 
stellar mass to an average massive galaxy in our sample from 
merging systems is highest around M* = 10 10 ' 8 Mq and de- 
clines at lower and higher masses. The gas mass accreted 
however is quite low for the highest mass galaxies, but in- 
creases at lower masses and stays relatively constant at val- 
ues M* < 10 10 - 5 M . 

Using this, we calculate how much gas mass is brought 
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Figure 5. The relative amount of gaseous and stellar matter 
accreted from mergers of various masses between redshifts 1.5 < 
z < 3. The x-axis is the stellar mass of galaxies merging with our 
typical average galaxy from our sample with M* = 10 112 Mq . 
The y-axis shows the relative contribution to the total for galaxies 
at masses separated by 0.05 dex. The solid line shows the stellar 
mass contribution and the dashed line the gaseous contribution. 



into our galaxies due to these minor mergers. We calculate 
that the average stellar mass weighted gas mass fraction is 
f g ~ 1.03 down to M* = 10 9 ' 5 M . From this, we calculate 
M g , M /M»(0) = M,, M /M*(0) x f g = 0.55±0.23. 

We can then use the above results and equation (9) 
to calculate the amount of gas accreted as a fraction of the 
initial stellar mass. Putting these together, we conclude that 



M g , A /M»(0) = 0.83 ±0.37. 



(13) 



What this implies is that a roughly constant gas mass frac- 
tion seen in the observations of these distant galaxies (§3.1) 
reveals that on order the entire initial stellar mass of a mas- 
sive galaxy is added over time outside of mergers to form 
stars during 1.5 < z < 3. This equates to an absolute 
amount of gas accreted as (1.8 ± 0.8) x 10 Mq over 2.16 
Gyr, or at an average rate of gas accretion of 



dM B , A (t) 
dt 



(83 ± 36) M Q yr" 



(14) 



needed to produce the star formation we observe. In the 
next section we evaluate the uncertainties within this mea- 
surement and in the section after give the implications of 
this result. 



4.3 Uncertainties 

There are systematic issues that could affect these results, 
some of which we examine in this section. One of these 
is that the the star formation density-gas density relation 
(Schmidt-Kennicutt law) may evolve (e.g., Daddi et al. 



2010). Since starbursts are more efficient, the initial gas 
mass could in principle be lower than what we are mea- 
suring. If the normal Schmidt-Kennicutt law then applied 
for the lower redshift galaxies in our sample the net result 
would be an increase in the gas mass to stellar mass ratio 
and equation 5 no longer strictly holds. Likewise, it may be 
possible that we are measuring the gas masses incorrectly 
too low at our highest redshift point for some reason and 
that there is actually a decrease with time in the gas mass 
fraction. We examine both of these scenarios quantitatively 
by considering how our relations would change if we use a 
gas mass density relation of the form: 



M g (t) 
M*(t) 



Mg(0) 
M*(0) : 



(15) 



where k is the relative change in the ratio between the gas 
to stellar mass at some lower redshift than at the initial time 
at z ~ 3. In this case, the relative gas mass accretion can be 
written as: 



M g , A (t) 
M,(0) 



, M E (0) <ib>8t 
kx (2.64T0.26) x + - 



M g , M (i) 



M»(0) M„(0) M,(0) 



(16) 



When using this equation with the value of k = 0.5, where 
the gas mass fraction has dropped by a half, we find that the 
gas accretion mass fraction drops to M g , A /M„ = 0.78 ±0.37. 
If we take the limiting case where all of the gas is exhausted 
from z = 3 to z = 1.5, where n — 0, we would find a lower 
limit accreted gas mass fraction of M g;A /M„ = 0.56 ± 0.37. 
In general, if the gas mass fraction decreases with time, with 
perhaps an evolving form of the relation between gas mass 
density and star formation rate density, this would result in 
a slightly lowered derived gas mass accretion. For the more 
likely case of less efficient star forming systems over time, 
the net gas mass fraction would increase even more than our 
initial calculation, given that in this case k > 1. However, 
most of the constraint on this measurement comes from the 
fact that there is a high ongoing star formation rate during 
this epoch, and the gas to fuel this star formation is not 
initially present within the galaxy, nor carried in through 
mergers. 

If we directly consider the more efficient Schmidt- 
Kennicutt law, as proposed by Daddi et al. (2010) for star 
bursting galaxies to apply throughout our redshift range, we 
would then have a lower gas mass fraction at both the start 
and the end of our evolution. As above, we would still have to 
account for the star formation present within these systems, 
and this is the driving observation behind our calculation 
of a high gas mass accretion rate. However it is important 
to note that our galaxies are in the regime where they fol- 
low the standard law (Daddi et al. 2010), where the more 
efficient starburst are for sub-mm and ULIRG galaxies. 

One remaining issue is that our star formation rates 
maybe too high, however, others have found very similar star 
formation rates or higher, for similar galaxies (e.g., Daddi et 
al. 2007). The star formation rates we measure are consistent 
in the ultraviolet and using Herschel far-infrared imaging 
(Hilton et al. 2012), and if anything we are at the lower end 
of the various star formation measurements (Bauer et al. 
2011). 
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4.4 Implications 

We present in this paper an analysis of the amount of cold 
gas mass which is likely accreted into massive galaxies with 
stellar masses M* > 10 11 Mq at 1.5 < Z < 3. This amount 
of accreted gas is necessary to account for the observed star 
formation rate between redshifts 1.5 < z < 3, which cannot 
be accounted for by gas brought in through minor+major 
mergers (Bluck et al. 2012). Our overall result is that we find 
a basic accretion rate of M g ,A = (83 ± 36) Mq yr -1 of cold 
gas needed from the IGM to account for this star formation. 

Papovich et al. (2011) recently published a measure- 
ment of the total gas accretion rate at similar redshifts, 
which includes all methods that bring in gas. They were 
not able to distinguish between the gas brought in through 
mergers, and that brought in through gas accretion. Their 
results are similar to ours, finding a higher accretion rate, 
but that the total accretion is similar to what we find here, or 
slightly higher. We note that we are only able to probe merg- 
ers for massive galaxies down to M* = 10 9 ' 5 Mq , and it is 
possible that the gas accretion we find originates from galax- 
ies with stellar masses lower than this limit. Furthermore, 
we are not able to say whether the accreted gas originates 
in clumps into the galaxy or in a smooth continuous mode. 
Our result is an integrated average over the time-spanned 
by redshifts 1.5 < z < 3. 

There are several implications from these results. The 
first is that what we actually measure is the amount of cold 
gas accreted which is turned into stars. Because of the sig- 
nificant amount of outflows in galaxies, particularly those 
with star formation rates similar to ours, a significant frac- 
tion of all accreted gas is removed from the galaxy by these 
winds (e.g., Dekel et al. 2009; Dave et al. 2011). The re- 
sult of these outflows, is that the amount of gas actually 
accreted is the combination of the outflowing gas as well as 
the gas used in star formation. This in principle could easily 
double the amount of gas mass needed to be accreted from 
the IGM (e.g., Faucher-Giguere et al. 2011). The ultimate 
morphologies of these galaxies in the local universe is also 
likely to be disky in some way, and would produce galaxies 
without large bulges as is found in merger remnants. This 
result also explains how massive galaxies in the very local 
universe, within 8 Mpc, can have almost no bulges (e.g., Ko- 
rmendy et al. 2010), as gas accretion is the dominant process 
for their formation. Gas accretion is also one way to solve 
the G-dwarf problem of having too many metal rich stars in 
the solar neighbourhood (e.g., Larson et al. 1974), and may 
relate to the gas accretion we see in our own galaxy (Blitz 
et al. 1999). 

Overall our result of finding an accretion rate of gas of 
M g ,A = (83 ± 36) Mq yr" 1 for the most massive galaxies 
at 1.5 < z < 3 is roughly consistent with theoretical ob- 
servations which predict a similar amount of gas accretion 
(Murali et al. 2002; van den Bosch 2002; Keres et al. 2005; 
Dekel et al. 2009; McBride et al. 2009). Early simulations 
by Murali et al. (2002) predict a gas accretion rate of M gi A 
~ 40 Mq yr" 1 , while more recent work suggests higher rates 
of Mg, A ~ 100 M Q yr" 1 (e.g., Dekel et al. 2009a,b; Faucher- 
Giguere et al. 2011). Our results are in general agreement 
with these models, although we are only examining a narrow 
redshift and stellar mass interval within this paper. 

Our results also reveal that amongst the material put 



into place in massive galaxies, the gas accretion accounts 
for 61±21% of the stellar matter added to galaxies from 
1.5 < z < 3. Mergers account for the remainder of the 39% 
of the mass added to galaxies over this epoch, with 1/3 of 
this minor mergers and 2/3 of this major mergers (Bluck et 
al. 2012). Overall this implies that gas accretion into massive 
galaxies at early epochs is a major formation method, and 
dominates over mergers as a formation process. 

Our results from this work, and previous papers examin- 
ing the issue of stellar mass formation modes (e.g., Conselice 
2006; Bluck et al. 2012) is similar to what is predicted in sev- 
eral studies using N-body simulations. This includes Genel 
et al. (2010), who find through simulations that about 60% 
of the halo mass in galaxies is put into place through merg- 
ing, with at least 40% of the mass coming from accretion. We 
also find about a 50-50 ratio of these two formation modes. 
Furthermore, using hydrodynamical models, Murali et al. 
(2002) find that for galaxies brighter than about a fourth 
of L* , the characteristic luminosity, gas accretion dominates 
over merging as a mechanism, whereby gas accretion, ac- 
counting for ~ 75% of the stellar mass build up at z ~ 2. 
While others, e.g., Stewart et al. (2008) find that over half of 
the mass in galaxies is formed in mergers, while others such 
as Angulo & White (2010) find that nearly all the mass in 
galaxies is put into place via mergers. These last few suggest 
that mergers are the dominant process, while we find that 
gas accretion is more important. However, our investigation 
is only for the most massive galaxies, and this may not be 
representative of typical galaxies at these redshifts. Future 
studies will have to probe deeper to obtain the merger his- 
tory for lower mass systems to carry out similar calculations 
as these and probe this evolution in different environments. 



5 SUMMARY 

We present in this paper a study of the cold gas mass densi- 
ties for massive galaxies with M, > 10 11 Mq at redshifts of 
1.5 < z < 3. While we do not directly detect the accretion of 
gas within our galaxies, we are able to make a circumstantial 
finding for its existence. Our conclusions are as follows: 

I. Utilizing our measured star formation rates and 
galaxy sizes we find a roughly constant cold gas mass to stel- 
lar mass fraction for this sample across the redshift range of 
1.5 < z < 3. 

II. We utilize the star forming and merging properties 
of these galaxies from previous work in Bauer et al. (2011) 
and Bluck et al. (2012) to measure the mass budget of our 
sample of massive galaxies, finding that A/ g .A = (83 ± 36) 
Mq yr -1 of gas is needed to sustain the star formation rate 
outside of gas brought in via mergers. 

III. We derive based on these values that cold gas ac- 
cretion from the intergalactic medium, or alternatively very 
minor galaxy mergers with mass ratios lower than 1/100, ac- 
counts for 61±21% of the baryonic matter added to galaxies 
from 1.5 < 2 < 3. 

This amount of gas mass added is larger than the 
amount of material added due to the merger process (both 
minor and major) (e.g., Conselice 2006; Bluck et al. 2012) 
and is largely in agreement with models which predict on 
the order of 100-200 M yr" 1 added from cold gas accre- 
tion (e.g., Dekel et al. 2009). Gas accretion is therefore the 



© 2002 RAS, MNRAS 000, \T§9\ 



Observational Evidence for Gas Accretion 9 



major method for building up the stellar masses of massive 
galaxies between redshifts 1.5 < z < 3. The remainder of 
the formation is due to merging. 

Further work with e.g., the CANDELS survey will allow 
us to carry out this measurement for lower mass galaxies 
where the mode of formation could be significantly different 
than the more massive systems (e.g., Dekel et al. 2009). Also, 
we are examining in this paper massive systems in the last 
throws of their formation at z < 3. Investigating the ratio 
of formation due to mergers and gas accretion at z > 3 will 
reveal how the first epochs of these massive galaxies were 
formed. This however will require observations from JWST 
and the ELTs. 

We thank the GNS team, particularly Fernando 
Buitrago and Amanda Bauer, for their contribu- 
tions to this survey and the previous published work 
utilised here. We thank the referee for a report that 
improved this paper significantly. The data and cat- 
alogs as used in the GNS survey are online at: 
http: / /www. nottingham.ac.uk/astronomy/gns/ 
The GNS is financially supported by STFC and the Lever- 
hulme Trust. Support was also provided by NASA/STScI 
grant HST-GO11082. 
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